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Abstract 

A convection  instability  characteristic  of  plasmas 
in  an  inhomogeneous  azimuthal  magnetic  field  is  treated  in 
the  linear  stage  and  in  nonlinear  saturation.  The  analysis 
is  done  in  such  a way  that  collisional  and  collisionless 
limits  can  be  taken,  and  these  limits  are  displayed  along 
with  the  more  general  Intermediate  result.  The  instability, 
known  previously  in  the  literature  in  its  collision-dominated 
form,^  is  shown  to  be  a "flute"  instability  with  collisional 
modifications  to  the  growth  rate.  The  nonlinear  saturation  is 
analyzed  by  examining  a finite  amplitude  restoring-force  term 
in  the  differential  equation  that  describes  the  instability. 
This  term  is  due  to  the  fact  that  the  Instability  convects 
plasma  into  striatlons  of  the  plasma  column  surface, 
modifying  the  density  gradient  until  the  restoring  forces 
balance  the  pressure  gradient  driving  force.  The  effects  of 
finite  ion  gyroradius  are  displayed,  and  applications  of  this 
study  to  convection  cells  in  a thermal  plasma  and  to  exploding 
wire  plasmas  are  discussed. 
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I.  INTRODUCTION 

The  conventional  approach  to  evaluating  the  stability  of  a 
plasma  configuration  makes  use  of  the  so-called  dispersion  relation. 

In  this  approach  one  starts  with  a number  of  first-order  differential 
equations,  with  time  as  the  independent  variable;  these  equations 
are  essentially  conservation  relations  which  describe  the  time  rate- 
of-change  of  physical  quantities  such  as  momentum,  energy,  density, 
etc.  Points  of  equilibrium  of  the  system  are  determined  by  setting 
the  time  derivatives  to  zero.  The  stability  of  an  equilibrium  is 
then  evaluated  by  applying  a small  perturbation  to  the  system  and 
setting  up  "variational"  equations  which  describe  the  temporal 
evolution  of  these  perturbations.  Typically  one  linearizes  these 
variational  equations  and  seeks  solutions  of  the  form  f (x)exp(iwt) , 

*\j 

where  the  frequency,  go,  is  a complex  quantity.  This  procedure  then 
usually  leads  to  the  formation  of  a determinant  which,  when  expanded, 
forms  an  algebraic  equation,  which  is  the  dispersion  relation.  The 
dispersion  relation  is  then  solved  for  go,  with  stability  depending  on 
the  imaginary  parts  of  the  various  roots.  While  it  is  useful  for 
determining  the  onset  of  an  instability  and  the  initial  growth  rate, 
the  dispersion  relation  is  not  always  convenient  for  analysis  of  the 
temporal  evolution  of  a perturbation  beyond  the  linear  stage. 

The  problem  of  evaluating  the  stability  of  a plasma  configura- 
tion can  also  be  treated  by  what  might  be  called  the  "different ial- 
-equation"  approach.  In  this  case,  one  again  starts  with  the  required 
first-order  differential  equations,  but  now  these  equations  are  combined 
to  form  a differential  equation  of  higher  order.  In  many  problems  of 
physical  interest  it  is  found  that  the  conservation  relations  can  be 
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a thermal  plasma  device  is  discussed  in  Sec.  VI  and  the  results 
of  the  work  are  discussed  in  Sec.  VII. 

II.  LINEAR  STABILITY  ANALYSIS 


Consider  a plasma-filled  gap  between  two  coaxial  cylinders  of 

radius  R and  R+d,  with  d<<R  (Fig.  1).  An  azimuthal  magnetic  field, 

which  falls  off  as  1/R,  exists  in  the  gap;  such  a field  might  be 

produced,  e.g.,  by  passing  a current  along  a conductor  located  at  the 

axis  of  the  cylinders.  The  surface  of  the  inner  cylinder  acts  as  a 

plasma  source,  for  example,  by  thermal  ionization'*'  while  the  outer 

cylinder  is  assumed  to  be  a sink  for  the  plasma.  It  will  be  shown 

that  this  system  is  unstable  since  the  density  gradient  is  parallel  to 

the  magnetic  field  gradient.  The  general  problem  is  treated  first 

and  the  collisional  and  collisionless  limits  are  then  compared  with 
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results  obtained  elsewhere. 

The  two-fluid  model  is  used  to  describe  the  plasma  and  the  iso- 
thermal equation  of  state  is  used  to  close  the  set  of  equations. 

The  quasineutrality  approximation  is  also  employed;  the  equation  of 
continuity  is 

u) 


where  the  subscript  ot=i,e  denotes  ions  and  electrons  respectively  and 
the  momentum  equation  is 


V-n 

3t  m c / m n t * 


\ Ta^n  la 


(2) 


where  and  r * are  the  collision  frequencies  of  neutrals  with  the 
ions  and  electrons  respectively.  Note  that  the  convective  derivative 
can  be  ignored  if  the  magnetic  field  is  sufficiently  large. 


expressed  by  algebraic  relations  rather  than  differential  equations, 
because  certain  physical  quantities  remain  constant  in  time.  Thus  the 
"order"  of  the  system,  which  is  equal  to  the  number  of  first-order 
differential  equations  required  for  its  description,  can  be  reduced 
significantly.  In  many  cases  the  overall  equation  reduces  to  a second- 
-order  differential  equation  with  time  as  the  independent  variable. 

It  is  then  possible  to  make  use  of  the  extensive  literature  that  is 
available  on  the  properties  of  these  equations  and  their  solutions. 

It  is  also  possible  to  identify  the  terms  which  play  the  role  of  restoring 
forces,  dissipation  etc.  and  to  determine  factors  that  enter  into 
linear  and  nonlinear  behavior;  these  capabilities  facilitate  the 
physical  interpretation  of  the  results.  The  role  of  different  time 
scales  also  becomes  clear,  so  that  standard  methods  of  singular 
perturbation  theory  can  be  used  in  obtaining  nonlinear  solutions 
in  which  such  time  scales  play  a role.  Finally,  this  approach  makes 
it  easier  to  evaluate  various  different  techniques  that  might  be 
employed  when  one  is  interested  in  stabilization  of  undesired 
instabilities  by  external  means. 

In  the  present  paper  the  differential-equation  approach  described 
above  is  used  to  analyze  an  instability  characteristic  of  a 

plasma  in  an  inhomogeneous  magnetic  field.  We  show  that  this 
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instability,  which  has  been  treated  earlier,  is  related  to  the  flute 
instability  in  a fully  ionized  plasma.  In  Secs.  II  and  III  of  this 
paper  the  linear  stability  analysis  is  carried  out  and  interpreted. 

Finite  Larmor  radius  corrections  to  the  pressure  tensor  are  discussed 
in  Section  IV,  and  in  Section  V the  nonlinear  saturation  of  the 
instability  is  analyzed.  The  possible  relevance  of  the  theory  to 
exploding-wlre  plasma  experiments  and  to  convection-cell  experiments  in 


Since  d<<R,  a local  cartesian  coordinate  system  with  origin 
in  the  outer  cylinder  can  be  used  as  an  approximation  as  shown  in 
Fig.  1.  Considering  only  the  dependence  on  the  x coordinate,  one 
can  write  the  steady-state  solution  of  Eqs.  (1)  and  (2)  as 
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Here  it  is  assumed  that  the  change  in  magnetic  field, 


j-S^,  across  the  gap  is  small  and  that  n(0)*0.  The  flow  in  the  -x 
direction  (radially  outward)  across  the  magnetic  field  is  due  to 
ambipolar  diffusion,  and  the  motion  in  the  z direction  arises  from  a 
combination  of  the  Jg*]^  drift  and  the  diamagnetic  drift.  Although 
the  gradient  B drift  drives  the  instability  it  is  much  slower  than 
both  the  E*J^  and  the  diamagnetic  drifts.  The  behavior  near  x=0  in  these 
equations  is  due  to  the  presence  of  a boundary  layer  of  thickness  on 
the  order  of  the  ion  gyroradius,  and  does  not  present  any  difficulties 
in  the  analysis. 

At  this  point  it  is  normally  assumed  that  the  unstable  oscillations 
in  the  system  will  occur  on  a time  scale  much  slower  than  x^  or  xg,  so 


that  the  Inertial  term  in  Eq.  (2)  is  usually  ignored. 
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However , 


that  procedure  reduces  the  resulting  equation  for  the  perturbed 
density  from  a second-order  differential  equation  to  a first-order 
differential  equation,  thereby  losing  one  of  the  roots  of  the  charac- 
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teristic  equation.  More  important,  keeping  the  inertial  term  leads 
to  corrections  to  the  growth  rate  of  the  instability  and  corrections 
to  the  critical  magnetic  field  at  which  the  instability  first  appears. 
The  tables  discussed  at  the  end  of  this  section  show  these  corrections 
and  the  ranges  in  which  they  are  important. 

Retaining  the  inertial  term,  we  take  the  cross  product  of 
Eq.  (2)  with  | and  solve  the  resultant  equation  for 

Substituting  this  expression  into  Eq.  (2)  yields  an  equation  for 
(the  symbol  "a"  means  that  a component  is  perpendicular  to  . 


2 

If  the  magnetic  field  is  sufficiently  large  then  (l/x^+S/St)  can  be 

2 

ignored  compared  to  in  Eq.  (7),  but  all  other  terms  are  retained. 

If  it  is  assumed  that  3/3y=0  (i.e.,  that  the  perturbations  are  azimu- 

thally  symmetric),  V can  be  set  equal  to  zero,  so  that  V =V 

’ay  n 'val 

When  the  inertial  terms  are  ignored  at  the  onset  of  the  calculations, 
the  time  derivatives  on  the  right  hand  side  of  Eq.  (7)  do  not  appear. 

This  procedure  is  equivalent  to  assuming  that  u<<l/r^,  whereas  Eq.  (7) 
extends  the  frequency  range  to  ioJvI/t^.  Thus  we  have 
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and  this  relation  holds  for 
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If  one  considers  electrostatic  perturbations  which  depend  only 
on  x and  z,  then  Eq.  (8)  can  be  substituted  into  the  continuity 
equation  (1),  and  the  resulting  equation  can  be  linearized.  The  next 
to  last  term  in  Eq.  (8)  is  the  only  term  which  does  not  contribute 
to  the  continuity  equation  when  is  considered  a constant  (since 
V.(JjxVn-)=0).  Thus  the  leading  order  contribution  of  this  term 
to  the  linearized  continuity  equation  results  from  the  presence  of  a 
magnetic  field  gradient.  The  effect  of  the  magnetic  field  gradient 
is  introduced  through  this  term  which  is  approximated  by 

cTa  _ /M”\  cT«  to' 


/P*vn  \ ci  . , 

a / ^ _ a 3n 

^ * \ = " eRB„  3z 


where  3/3x[Bq(^) ]=-1/RBq.  Here  the  assumption  d<<R  has  been  used. 
Again  if  Bq  were  truly  constant  then  this  term  would  vanish.  Note 
that  a prime  signifies  a first-order  perturbed  quantity. 

4 

Following  the  method  used  by  Kadomtsev  and  Nedospasov,  and 
Simon^  the  x dependence  of  the  perturbed  quantities  is  expressed  by  a 
Fourier  sine  series  in  which  only  the  first  term  is  retained;  the 
linearized  continuity  equation  is  then  averaged  over  x from  0 to  d 
weighted  by  sin(irx/d).  This  procedure  yields  two  coupled  first  order 
differential  equations  involving  the  perturbed  density,  n' , and 
the  perturbed  electric  potential,  <(» * . These  equations  can  be 
combined  to  obtain  a single  second  order  differential  equation  in 
the  perturbed  density: 

d-p-— t-  + A + A n'(t)  = 0 (11) 

dt2  1 dt  2 

where  now  n'  is  only  a function  of  time  and  where 
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Here  k is  the  axial  wavenumber  and  k=tt  /d  +k  is  the  total  wavenumber. 
The  numerical  factors  and  which  are  obtained  by  averaging  over 
the  x coordinate,  are  given  in  Appendix  I. 

Since  these  are  very  complicated  expressions,  the  coefficients 
will  be  simplified  by  specializing  to  the  case  T^=Te=T  and  by  using 


the  Spitzer  result 


(fir)i  (x/m)i  ^m£ 
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<<  1 . 


For  simplicity  the  remainder  of  this  section  will  pertain  to  this 
special  case.  Under  these  conditions  A^,  A£  and  8 reduce  to  (n  is 
a numerical  factor  and  is  defined  in  Appendix  I) 
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where  the  last  approximate  equality  in  Eq.  (18)  holds  as  long  as  K r^<l 
(r^  is  the  ion  Larmor  radius) . In  Section  IV  we  will  show  that  this 
is  a reasonable  assumption.  Also  note  that  E^,  which  is  related  to  the 
ambipolar  electric  field  through  Eq.  (4),  reduces  to  T/e  in  this 
special  case.  In  addition  the  axial  drift  velocities  given  by 
Eq.  (6)  reduce  to 


V,  « 0 
iz 


ez  eBx 


For  the  ions  the  diamagnetic  drift  essentially  cancels  the 
drift  whereas  for  the  electrons  the  two  drifts  add. 

Since  and  A^  are  complex,  the  solution  of  Eq.  (11)  will 
exhibit  oscillatory  behavior  in  the  laboratory  frame  of  reference. 

This  oscillatory  behavior  results  from  the  drift  motion  of  the  electron 
and  ion  fluids.  The  complex  characteristic  roots  of  differential 
equation  (11)  with  coefficients  A^  and  given  by  Eqs.  (16)  and  (17) 
are  easily  found  to  be 
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If  the  critical  magnetic  field,  B^,  is  defined  as 
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2 2 

then  for  B /B  >1  the  two  roots  given  in  Eq.  (19)  simplify  to 
c ^ 
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where 
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From  Eq.  (19)  it  is  clear  that  the  part  of  which  determines 


the  oscillation  frequency,  is  related  to  the  average  drift  velocity  of 
the  electron  fluid,  so  that  the  z dependence  of  the  perturbation  may 
be  written  as 


n1(x,z,t)  = n1(x,t)e 


ik(z-vDt) 


(24) 


with 


_ 2nTc 
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(25) 


The  second  mode,  u^,  can  be  associated  with  the  ions,  which  are 


stationary  in  the  laboratory  frame,  so  that  is  a purely  damped 
mode.  This  result  is  similar  to  the  case  of  the  universal  over- 
stability of  a resistive,  inhomogeneous  plasma  where  one  mode  is 


a purely  growing  mode  in  the  frame  drifting  with  the  electrons,  and 
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the  other  mode  drifting  with  the  ions  is  damped. 

If  one  transforms  equation  (11)  into  the  frame  drifting  with  the 


electron  fluid,  coefficients  and  become 
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Now  if  A ^ were  merely  equal  to  1/t^,  it  would  represent  a positive 


dissipation  term  and  could  not  lead  to  Instability.  On  the  other 

hand,  A2,  which  represents  a restoring  force  term,  can  lead  to  instability 
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if  it  becomes  negative  (i.e.,  if  B>Bc> . Ip  fact  Eq.  (11)  with 

real  and  positive  and  real  and  negative  will  be  recognized  as  the 

equation  for  an  inverted  pendulum  and  it  can  be  shown  that  the 

general  class  of  "convection"  instabilities  such  as  the  Benard 
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instability  can  be  described  by  an  equation  of  this  form.  ’ 

This  class  of  instabilities  is  characterized  by  purely  growing  "zero- 
frequency"  perturbations  which  arise  from  a negative  restoring  force. 
This  present  problem  is  very  similar  in  that  the  instability  is 
driven  by  a negative  restoring  force,  but  more  general  since 
is  complex. 

For  the  collision  dominated  case,  where  Uq«y<1/t^,  ( y is  the 
growth  rate)  and  in  the  collisionless  limit  (Wq.y^I/t^)  Eq.  (11) 
can  take  the  form  of  the  inverted  pendulum  equation  with  both  A^ 
and  A ^ real.  In  the  laboratory  frame  for  this  collision  dominated  case 
the  coefficients  become 
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Here  the  drift  moM.on  occurs  on  a time  scale  much  slower 


growth  of  the  instability.  For  the  collision  free  case, 
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drifting  with  velocity  the  coefficients  become 
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In  both  cases  the  general  results  derived  for  the  "convection"  type 
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instability  can  be  applied  directly.  * From  Eq.  (31)  it  is  clear 
that  the  critical  magnetic  field  must  be  redefined  (see  Table  III) 
for  the  collisionless  case. 

Working  again  in  the  laboratory  frame  Table  I summarizes  some 
of  the  results  of  Eq.  (19)  for  the  unstable  mode.  The  oscillation 
frequency,  id^,  and  the  growth  rate,  y,  are  listed  for  the  cases  with 
B>Bc  and  with  where  B must  be  greater  than  Bc  in  order  for 

instability  to  occur.  Note  that  even  in  the  laboratory  frame  the 
instability  is  essentially  a purely  growing  mode  in  the  very  large 
magnetic  field  limit  (w^O,  y>0) . The  growth  rate  of  the  instability 
is  plotted  as  a function  of  the  axial  wavenumber,  k,  in  Fig.  2,  for  the 
case  where  B>>Bc.  From  this  plot  we  see  that  the  growth  rate  is 
largest  when  k>ir/d.  Wavenumbers  much  larger  than  rr/d  should  be 
stabilized  by  finite  Larmor  radius  corrections  to  the  ion  pressure 
tensor.  This  point  is  discussed  in  detail  in  Sec.  IV. 

In  the  low  frequency  collisional  limit  where  l/x/»  |w |= j w^+iy | , 
the  results  of  Eq.  (19)  are  found  in  Table  II.  In  the  collisionless 
limit  the  results  of  Eq.  (19)  are  listed  in  Table  III. 

The  results  listed  for  the  low  frequency  collisional  limit 

(Table  II)  agree  with  those  found  by  Timofeev.^-  It  is  not  possible 

to  directly  establish  from  Timofeev's  results  which  collisionless 

instability,  if  any,  is  the  analog  of  the  collision;.l  instability  since 

he  assumes  l/x^» | w | = | Wg+iy | ; however,  since  the  expression  for  the 

growth  rate  in  the  collisionless  limit  (Table  III)  agrees  with  the 

expression  for  the  growth  of  this  type  of  flute  instability  found  by 
3 

Krall,  it  can  be  inferred  that  the  instability  described  in  Tables  I 
and  II  is  the  collisional  analog  of  the  collisionless  flute  instability. 


- .. _ 
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Note  that  when  B>>Bc,  the  growth  rate,  Yq,  for  the  collisionless  case 
(Table  III)  is  related  to  the  growth  rate,  yc>  for  the  collisional 
case  (Table  II)  in  the  usual  way: 


2 

Y0 


c v 


coll 


Vi 


Figure  3 shows  the  growth  rate  as  a function  of  the  collision 
frequency,  r \ f°r  the  case  k>>ir/d  and  B>>Bc- 

III.  INTERPRETATION  OF  INSTABILITY  CRITERION 

The  interpretation  of  the  instability  criterion  is  not  to  be 
found  simply  in  collisional  damping  of  the  mode,^  as  one  can  see 
by  taking  the  collisionless  limit  of  Eq.  (19).  This  gives  the 
condition  (Table  III): 


I±  < A 

,2  ( 22  V 

d ' it  n ' 


(33) 
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which  is  the  condition  for  finite  gyroradius  stabilization  of 

the  interchange  instability  with  driving  acceleration 

2 


g 


thi 

2n2R 


(34) 


and  typical  wavenumber  k^tr/d,  as  found  in  this  paper.  This  is  the 


dominant  criterion  until  collisions  become  frequent  enough  that 

1/4 


2Kd 


i ^ qkv 


thi 


(5) 


(35) 


i.e.,  until  the  mean  free  path  becomes  shorter  than  d.  When  this 

occurs,  the  instability  criterion  is  (Table  I) 

1/2 


(fit)2  > 


(36) 


'I 


.1 


for  k'vir/d,  rather  than  the  condition  given  in  Eq.  (33).  In  summary, 


one  must  have 


i.e. , 


d/R  > Max 


d/R  > Max 


2 d2 

2*2  ( 

m 1/2 
m x 

— ) 

(«T)2  ' 

m / 

/ 2 2 

| 

If 

R2 

1 2 

m 1/2 
(A 

2 IT 

(flT)2 

\ m^/ 

(37) 


for  instability.  Steeper  density  gradients  are  stabilized  by 

finite  mean  free  path  and/or  finite  gyroradius. 

When  the  magnetic  field  is  self-generated,  as  in  a current- 

carrying  plasma  pinch,  one  should  compare  Eq.  (37)  with  the 

12 

instability  condition 


-d(fcnP) /d(£nr)  > 2y 


(38) 


(■^adiabatic  exponent  and  should  not  be  confused  with  the  growth  rate) , 


i.e. , 


d/R  < (2y)  A , 


(39) 


obtained  from  the  energy  principle 

6P5U  + ^ (6U) 2 < 0 (40) 

for  instability.  (Here  U = jdfc/B  and  P is  the  plasma  pressure.) 

This  latter  condition  sets  a maximum  value  for  unstable  d's,  while 
Eq.  (37)  sets  a minimum  value. 
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IV.  FLR  CORRECTIONS  TO  THE  PRESSURE  TENSOR 


The  addition  of  finite  Larmor  radius  (FLR)  corrections  to  the 
ion  pressure  tensor  renders  the  mathematics  very  cumbersome  and  was  not 
attempted.  However,  insight  can  be  obtained  by  comparing  the  magnitude 
of  such  correction  terms  with  the  magnitude  of  a typical  term  in 
the  linear  analysis  of  Section  II.  If  the  FLR  correction 
terms  are  found  to  be  significant  in  this  comparison,  one  would 
expect  FLR  stabilization  of  the  instability. 

The  magnitude  of  a typical  FLR  correction  term  to  Eq.  (2)  is  given 
by  r2ft^/4  S2V,  where  S2«MAX(k,  ir2/d2).^  Thus  in  the  colllsional 


regime  FLR  corrections  are  important  if 


(r^is2V/4) 


(flT>i  » 1 


(41) 


Similarly  in  the  collisionless  regime  the  condition  is 


(r^±S2V/4) 
(V/  | co  | ) 


2 2n 

s rA  > 

'V, 


4 M 


1 


(42) 


In  Figure  4 the  critical  magnetic  field  for  instability  is  plotted 
as  a function  of  the  axial  wavenumber,  k,  for  the  colllsional  regime. 
Here  (ftx)^1^*  is  given  by 

9 m 1/2  1/2 

- ¥ (£(<)  ) 

Th’  dashed  line  separates  the  regions  where  the  FLR  corrections  are 

large  and  small.  A similar  plot  for  the  collisionless  regime  is 

— crit 

shown  in  Figure  5.  In  this  case  * is  given  by 

<I!iVcrlt  ■ '?  <44> 


where  t Is  a convenient  time  scale  defined  by 
N 
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i(V-JL) 

*>'  T nV 


1/2 


(45) 


Unlike  the  collisionless  case,  collislonal  stabilization  dominates  for 
k<ir/d  in  the  collislonal  regime,  whereas  FLR  stabilization  dominates 
for  k>ir/d  in  both  the  collislonal  and  collisionless  regimes. 

The  largest  axial  wavenumber,  k^,  which  is  unstable  is  determined 
by  evaluating  Eq.  (41)  (or  Eq.  (42))  at  the  critical  magnetic  field 
given  by  Eq.  (43)  (or  Eq.  (44)).  This  marginally  unstable  k^  value 
corresponds  to  the  point  of  Intersection  of  the  solid  and  dashed  curves 
in  Figure  4 (or  5) , and  the  shaded  area  denotes  the  region  of 
instability.  The  FLR  corrections  can  actually  completely  quench 
the  instability  if 


£ 

d 


32  r~ 


Rd 


1/2  1/6 


(?)  ) 


(46) 


for  the  collislonal  case  or 


TT 

d > 


<Rr*>1/3 


(47) 


for  the  collisionless  case.  From  this  it  is  clear  that  Figure  2 

should  be  modified  particularly  for  large  values  of  k where  FLR 

stabilization  is  most  Important.  In  fact  the  curve  in  Figure  2 

should  fall  off  such  that  the  value  of  k corresponding  to  the  maximum 

growth  rate,  k , should  fall  between  -7  and  k_ . 

max  a L 


£ 

d 


< 


k < kT 
max  L 


IV.  NONLINEAR  SATURATION 


(48) 


The  understanding  and  identification  of  nonlinear  mechanisms 
capable  of  limiting  the  growth  of  the  instability  described  above 
is  facilitated  by  using  the  differential-equation  approach.  It 


has  been  argued  in  connection  with  Eq.  (11)  that  the  Instability 
arises  from  a negative  restoring  force  due  to  the  presence  of  a 
density  gradient  parallel  to  the  gradient  of  the  magnetic  field. 

This  feature  suggests  that  saturation  of  the  instability  can  be 
obtained  from  a mechanism  that  provides  a nonlinear  term  corres- 
ponding to  a positive  restoring  force  large  enough  to  overcome  the 
negative  restoring  force  in  Eq.  (11).  It  can  be  shown  that  the  simplest 

g 

term  that  meets  this  specification  leads  to  an  equation  of  the  form 

^ + A2y  + A3y3  = 0 (49) 

where  A3>0,  A^>0  and  A2<0.  Evidently  this  equation  is  satisfied 

g 

by  the  saturation  amplitude 


Since  the  coefficients  in  Eq.  (49)  are  all  real  this  quasilinear 
analysis  will  strictly  only  apply  to  the  collision  dominated  case 
described  in  Eqs.  (28)  and  (29),  and  to  the  collisionless  case  described 
in  Eqs.  (30)  and  (31);  it  is  expected  however  that  the  more  general 
problem  with  complex  coefficients  would  have  similar  results. 

In  considering  the  physics  of  the  problem  it  is  noted  that 
the  gradient  of  the  magnetic  field  can  not  be  affected  by  the 
instability  because  of  the  electrostatic  nature  of  the  Instability. 
Therefore  it  is  anticipated  that  the  saturation  effect  will  arise  from 
a modification  of  the  density  gradient  caused  by  the  instability  itself. 
Once  the  nonlinear  component  of  the  restoring  force  has  been  identified 
the  saturation  amplitude  can  be  calculated  from  Eq.  (50)  and  its  effect 


on  the  linear  density  gradient  can  be  determined. 
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In  order  to  derive  an  expression  for  the  nonlinear  term  in  Eq.  (49), 
Eq.  (8)  is  used  in  the  continuity  equation,  which  is  not  linearized 
in  this  analysis.  Again  an  expression  similar  to  Eq.  (10)  is  used 
to  describe  the  gradient  in  the  magnetic  field.  If  we  first  write 


where 


n(x,z,t)  = nNL(x)+n' (x,z,t) 


<j>(x,z,t)  = 4>NL(x)4-^  (x,z,t) 


n^L(x)  * n°(x)  + <5n°(x) 


4>^(x)  =>  <*>°(x)  + 6(j>° (x) 


(51) 


(52) 


the  continuity  equations  (one  for  each  species)  can  be  integrated  over 
9 

the  z coordinate.  Here  n and  4> ' ate  periodic  in  z in  the  frame 
drifting  with  the  electron  fluid  and  represent  the  perturbation  due 
to  the  unstable  mode.  A time  average  is  not  appropriate  since  the  pertur- 
bation is  purely  growing  in  the  drift  frame.'1  Note  that  6n^  and  6<t>®  are 
the  average  nonlinear  modifications  to  the  equilibrium  density  and 
electric  potential  respectively. 

After  integration  over  z the  steady  state  equations  are  (a*i,e) 


^7  • ^ |^7n°+en'b  * 7$’  - (—  + |r)  7<f>’ 

N'v.  ft  02  9t  'v  ^ ,V/T  ft  It  3t / <v< 

a n a ' a ' 


(53) 


0 «~4n 


(St)  " (ftx)  l*' 


q T 2 0 
+ -JLJL- _ 0 


z e(QT)a  3x2 


The  only  time  dependence  in  this  equation  arises  from  the  drift  motion 
since  in  steady  state  the  amplitude  has  saturated;  therefore  u)“u)Q+iY»WQ 
in  what  follows.  Since  6n^  and  64P  are  assumed  to  be  small  corrections 
to  n^  and  4?  respectively,  terms  of  the  form  6n^6^  are  ignored. 
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Eliminating  <5$^  and  solving  for  6n^  results  in  the  equation 


326n° 

3x2  “Te 


ir  <Z  ‘ Kwi*  ^lf^n0-w((nT)i+(nT)e)n,W 

i v n 


+ en,(VTi)^(|f)  > 

z 

The  following  expressions  are  used  for  n'  and  4>'  (the  subscript  s 
signifies  the  saturated  value  of  the  amplitude) 


n'(x,z,t)  - n'sin(kz-u).t)sin(— ) 
s 0 a 


<t>'(x,z,t)  « [ d> ’ sin(kz-a)  t)  + cos(kz-o).t)  ]sin(— ) 

s,s  u s,c  U a 

where  and  <f>'  allow  for  the  difference  in  phase  between  n'  and 

S | S S | c 

t ' . Solving  Eg.  (54)  for  6n^  yields 

e[((lT)  +<(!,)  1 

5n  — fir- — 

e i 

(56) 

e(T  -t_)  un  . . «. 


i e , 0S  / ,2irxs  , \ 

(T>  ”s*s,c  ; 1 - V 


where  the  boundary  condition  6n  (d)=0  has  been  invoked. 

In  order  to  relate  n'  and  <f>'  , the  same  procedure,  used  to 

S S j c 

obtain  Eq.  (53),  is  repeated  with  the  exception  of  multiplying 
by  n'  before  integrating  over  z.  The  resulting  equation  is 

0 - (T.+Ijln;  sin  (-j-)  - ( -j ) <— >°X,csl"  <T> 


('n',i)”^i.c’,8l"2(T)  +(-ix)(-§-)  nX,c8l"(T,c»s(f> 


As  in  the  linear  analysis  of  Section  II,  the  appropriate  procedure  at 

this  point  is  to  average  over  x since  the  functional  dependence 

on  x is  only  approximate.  Integrating  Eq.  (57)  over  x and  solving  for 
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s,c 


obtains 


n <(> 
s s,c 


d(T  +T  )n'2 
e i s 


(nT)e+(0T)1  V^e'V 


eNk 


4k 


so  that  fin  becomes 


n d2n'2[(Ox)  +(0x).]sin(^) 
6n°  = — e i d 


8irN 


(fix)  +(fix),  und 

e i . U / \ 

-2 + 4T  <V'i>, 


(58) 


(59) 


V,2d(T1-Te)  (cos(^)~l) 


8kN 


(-(fix)  +(fix)  o>nd 

e i , u x 

"7 + 4F 


If  one  considers  the  case  where  k=ir/d  and  Te=:Ti,  then  Eq.  (59) 
simplifies  to 


0 7f(ns)  f2irxx 
6n  = ~ 4N sin(~> 


(60) 


Here  the  assumption  that  used  in  Eq.  (7),  is  also  invoked, 

although  here  may  be  shifted  from  the  linear  value.  This  modification 
of  the  initial  density  gradient  due  to  the  presence  of  the  instability 
is  shown  in  Fig.  6.  Note  that  the  modification  is  such  as  to  decrease 
the  density  gradient  in  the  interior  of  the  plasma  and  thus  shut  off 
the  instability. 

If  this  nonlinear  mechanism  is  incorporated  into  the  previous 
analysis,  the  significant  modification  arises  from  the  additional 
nonlinear  term  in  the  restoring  force,  resulting  in  an  equation  of  the 
form  of  Eq.  (49).  Using  Eqs.  (13),  (14),  (50)  and  (60),  one  finds  that 
the  saturated  amplitude  of  the  density  perturbation  is  given  by 
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, 2N  A 
1 ■ — (1 
s n \ 


d(fir)  (fir), 
e 1 


In  the  collisional  limit  Eq.  (61)  can  be  written  as 

/ B2  \1/2 


which  is  valid  as  long  as  B>Bc.  If  B becomes  too  large,  however,  the 
nonlinear  analysis  is  somewhat  suspect,  since  ng  was  assumed  to  be  small 
compared  to  N.  In  the  collisionless  limit  Eq.  (61)  simply  becomes 


Here  the  analysis  is  somewhat  suspect  for  all  B>Bc.  However  it  is 
clear  that  the  amplitude  can  grow  much  larger  in  the  collisionless 
case,  than  in  the  collisional  case. 

It  should  be  pointed  out  that  there  may  be  other  important  nonlinear 
processes  which  have  not  been  included  in  this  model. 


VI.  EXPERIMENTAL  APPLICATIONS 


The  plasma  configuration  used  as  a basis  for  the  present  instability 
analysis  can  be  realized  physically  in  an  experimental  arrangement 
that  exploits  the  "vapor-pressure"  mode  sometimes  used  in  the  Q-machine.1^ 
In  this  case  the  role  of  the  inner  cylinder  in  Fig.  1 is  played  by  a 
metal  cylinder  of  a refractory  metal  like  tungsten,  which  is  heated  to 
an  electron-emitting  temperature  (2000°K).  This  cylinder  is  used  for 
thermal  ionization  of  an  appropriate  vapor  and  acts  as  the  plasma  source. 
The  outer  cylinder  in  Fig.  1 is  a heat-resistant  dielectric  cylinder. 

This  cylinder  acts  as  a plasma  sink  since  the  outwardly  drifting 
plasma  which  reaches  it  is  lost  by  recombination  at  the  surface. 


The  annular  space  is  flooded  with  Cg  or  K vapor.  An  azimuthal  magnetic 
field  which  falls  off  as  1/r  can  be  produced  in  a straightforward 


22 


way  by  passing  a current  through  a single  straight  conductor  along 

the  cylinder  axis,  and  calculations  show  that  the  required  fields 
2 3 

(10  -10  g)  can  be  produced  conveniently  by  conventional  power  sources. 
The  time  and  space  dependence  of  relevant  plasma  quantities  such  as 
the  density  can  be  measured  by  wall  probes  installed  so  as  to  be  flush 
with  the  surface  of  the  outer  cylinder,  thus  minimizing  plasma 
perturbations.  An  experimental  arrangement  of  this  kind  is  of  interest 
in  that  it  would  be  possible  to  observe  the  time-space  development 
and  structure  of  the  convection  cells  characteristic  of  the  instability. 
This  arrangement  would  also  be  useful  for  the  evaluation  of  various 
stabilization  techniques. 

Another  possible  application  of  the  formalism  for  this  instability 

is  to  the  plasma  formed  by  exploding-wire  discharges.  In  the  corona 

of  such  a plasma,  one  can  have  parallel,  inwardly  directed  gradients 

of  plasma  density  and  cf  the  azimuthal  self  magnetic  field.  Depending 

on  the  local  temperature,  the  plasma  at  these  radii  may  be  collisional 

or  not.  At  typical  temperatures  of  order  200  eV  and  densities  of  order 
19 

10  , collisionless  dynamics  are  appropriate.  At  this  temperature, 

deuterium  has  gyroradius 


r^CcmJ^/BCkG)  . (64) 

The  self  field  can  be  in  the  megagauss  range,  so  that  the  condition 

r 2 

B > Bc,  i.e.,  (-^)  « | (65) 

can  be  satisfied  for  low-atomic  number  plasma  (e.g.,  deuterium)  even 

when  d is  quite  small  (e.g.,  density  falloff  over  a distance  of 

-2 

d'v-lO  cm.)  The  formation  of  beads,  or  ripples  in  the  radius,  of 
such  plasmas  is  a well-known  phenomenon  and  is  much  more  violent  for 


14 

low  atomic  numbers.  If  this  instability  is  the  cause,  a typical 
wavelength  on  the  order  lcwr/d  should  be  expected,  since  finite 
Larmor  radius  effects  reduce  the  growth  rates  at  higher  k values, 
d is  estimated  at  0.05  cm  and  (as  defined  in  Sec.  IV)  is  ^250  cm  \ 
so  that  Eq.  (48)  predicts  that  the  wavenumbers  corresponding  to 
maximum  growth  is  in  the  range  60  cm  cm  The  observed 

periodicities  range  from  k^60  cm  ^ to  kvl50  cm  ^ . The  conditions 
for  the  ordinary  sausage  instability  in  this  current  carrying 
plasma  are  probably  also  satisfied  in  such  plasmas;'*''’  for  200  eV 
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electron  temperatures,  Bge^^10  gauss  and  n^lO  , the  sausage  and 
flute  growth  rates  are  comparable,  the  sausage  growth  scaling  as  B 


-1/2 


and  the  flute  as  T 


1/2 


DISCUSSION 


In  this  paper  a nonlinear  second-order  differential  equation 
has  been  used  to  analyze  the  linear  onset  and  nonlinear  saturation  of 
a convection  instability  characteristic  of  an  inhomogeneous  plasma  in 
an  inhomogeneous  magnetic  field.  This  differential-equation  approach 
has  been  used  to  extend  the  frequency  range  of  the  usual  dispersion- 
relation  analysis  to  the  collisionless  limit  and  the  linear  stability 
analysis  shows  the  instability  to  be  the  collisional  analog  of  the 
flute  instability,  both  being  driven  by  a density  gradient  parallel 
to  a magnetic  field  gradient.  A critical  magnetic  field,  Bc,  is  predicted, 
with  the  plasma  becoming  unstable  when  B>Bc*  When  B<Bc  the  plasma  is 
stabilized  by  finite  gyroradius  effects  in  the  collisionless  case  and 
by  the  effect  of  the  mean-free-path  in  the  collisional  case.  Finite 
gyroradius  corrections  to  the  pressure  tensor  are  also  found  to  be 
important  for  large  values  of  the  wave  vector  perpendicular  to  the 


24 


magnetic  field. 

The  form  of  the  nonlinear  differential  equation  indicates  that 
the  present  instability  belongs  to  a class  of  instabilities 
like  the  Benard  instability  of  a thermally  unstable  fluid  layer 
and  oth^r  instabilities  that  can  be  described  by  the  equation  for  an 
inverted  pendulum.  These  are  "symmetry-breaking"  zero-frequency 
instabilities  in  which  a medium  that  is  spatially  uniform  in  the 
stable  state  exhibits  a periodic  spatial  structure  consisting  of 
convection  cells  in  the  unstable  state.  The  ultimate  size  of  such 
convection  cells  is  determined  by  a spatial  average  (over  a cell) 
of  the  dynamic  energy  balance,  which  yields  the  saturation  amplitude 
of  the  instability. 

The  theory  of  this  instability  has  been  applied  as  a possible 
additional  mechanism  to  the  sausage  instability  to  explain  the  "beading" 
effect  observed  in  exploding-wire  plasmas,  and  a thermal-ionization 
experimental  arrangement  has  been  described  for  the  experimental 
observation  of  the  plasma  convection  cells. 
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APPENDIX  I 

In  performing  the  average  over  the  x dependence  to  arrive  at 
Eq.  (11)  a number  of  numerical  factors  are  encountered.  They  are 


iABLE  I - Growth  rate,  y,  frequency,  u>_,  and  critical  magnetic 
field,  Bc>  for  instability.  K*(k^-Hr2/(j2jl/2  tj,e  total 
wavenumber,  r^  is  the  ion  Larmor  radius  and  v^j  is  the  ion 
thermal  velocity,  ri  is  defined  in  Appendix  I. 
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TABLE  II  - Growth  rate,  y,  frequency,  u>q,  and  critical  magnetic 
field,  Bc,  for  instability  in  the  collision  dominated  limit. 
K»(k^-Hr^/d^)^'  ^ tjje  total  wavenumber,  r^  is  the  ion  Larmor 
radius  and  vthi  is  the  ion  thermal  velocity,  n is  defined  in 
Appendix  I. 
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TABLE  III  - Growth  rate,  y»  frequency,  wq,  and  critical  magnetic 
field,  Bc,  for  Instability  in  the  collisionless  limit.  K« 
(k2+ir^(i2)l/2  is  tjje  totai  wavenumber,  r^  is  the  ion  Larmor 
radius  and  vt«  . is  the  ion  thermal  velocity,  n is  defined 
in  Appendix  I. 


FIGURE  CAPTIONS 


Fig.  1 


Fig.  2 


Fig.  3 


Fig.  4 


Fig.  5 


Fig.  6 


Plasma  filled  gap  between  two  coaxial  cylinders  of  radius 
R and  R+d  (R>>d).  The  azimuthal  magnetic  field  and  a local 
Cartesian  coordinate  system  are  shown. 


Plot  of  growth  rate  as  a function  of  the  axial  wavenumber,  k. 

For  the  case  shown  here  B»B  and  T*T  =T;  also  yU4V  = - 

/2  c i e MAX 


_!.+ 1(-L  + ! §LX' 

2t4  2 V 2 dRm, / * 


Plot  of  growth  rate  as  a function  of  collision  frequency  with 

k>>7r/d.  For  the  case  shown  here  B>>Bn  and  T =T  ; also  y = 
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Plot  of  the  critical  magnetic  field  (solid  curve)  as  a 
function  of  the  axial  wavenumber,  k,  in  the  collisional 
regime.  The  dashed  curve  separates  the  regions  where  the 
FLR  corrections  are  large  and  small.  The  shaded  area  denotes 
the  region  of  instability  (Qt)*IN  = ] . (fiTi>MAX= 

/ j2  i 

~2  2 »and  a is  defined  in  Eq.  (33). 
n ri 

Plot  of  the  critical  magnetic  field  (solid  curve)  as  a function 
of  the  axial  wavenumber,  k,  in  the  collisionless  regime. 

The  dashed  curve  separates  the  regions  where  the  FLR  corrections 
are  large  and  small.  The  shaded  area  denotes  the  region  of 


instability,  is  defined  in  Eq.  (32)  and 
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Plot  of  the  density  as  a function  of  position  in  the  linear 
stage  (labeled  by  n^(x))  and  in  the  nonlinear  stage  (labeled 
by  n^(x)). 


